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The snake phenomenon has been observed in JET during the experiments 
{Sj ■ of pellet injection and consists in formation of persistent density perturbations 

(3jT), at rational-g surfaces. These structures persist over several sawtooth collapses 

and are difficult to explain as magnetic perturbations. Possibly related to this, 
there are indications that the tokamak plasma density has an anomalous radial 
pinch, much larger than that of the neoclassical origin. In this class of phe- 
nomena one should also include the persistent impurity accumulation in laser 
blow-off injected impurity, observed in experiments in TCV. There are several 
studies of the statistical properties of the correlations between the peaking 
\ factors for density, current density or pressure, with plasma parameters and 

these studies seem to support the idea of turbulent equipartition of the thero- 
modynamic invariants (Minardi). However we should note that these studies 
involve quantities expressed as global variables (like averages) and they can 
q ! hide other dependences not immediately obvious. 

We consider the possibility that the particle density behavior (and partic- 
ularly the snake phenomenon) can be connected with the existence of attract- 
ing solutions of certain nonlinear integrable equations. The statistical studies 
carried out on a large set of discharges with the purpose of testing the pre- 
7— i ! diction of the Turbulent Equipartition theory have suggested that the current 

density is given by the equation Aj + (A 2 /4) j = where A is a constant. 
However this may be valid on finite spatial patches and on every patch the 
equation can be considered an approximation of some more general equation, 
for which the derivation from first principles may be possible. One is the sinh- 
Poisson equation, A<p + (A 2 /4)sinh0 = for which the previous equation is 
the linearised form. The reason to consider the sm/i-Poisson equation comes 
from the existing proofs that this equation governs the asymptotic states of 
ideal fluids, or, more generally, of 2D systems that can be reduced to the dy- 
namics of point-like elements interacting by the potential which is the inverse 
o ,1 of the Laplacean operator (Jackiw and Pi, Spineanu and Vlad). This equa- 

j> | tion is however obtained when there are two kinds of elements (like positive 

and negative vorticity) and they are of equal number, n + = Then the 



X 



sm/i-Poisson equation is obtained as governing the states with maximum en- 



tropy of the discrete statistical system at negative temperature (Montgomery 
et al). Since the equation for the current density mentioned above is de- 
rived under the assumption of turbulent equipartition, the two descriptions 
may be related. However, the solutions for the unbalanced system of ele- 
ments, a = ^ 1, A(f> = (A 2 /8) [exp (0) / y/a — A/aexp (—0)] have been 
obtained numerically (Pointin and Lundgren) and have been shown to have 
higher entropy and higher stability than those of the sm/i-Poisson, precisely 
the characteristics we are seeking for. The limiting form of the unbalanced 
equation is the Liouville equation, A0 = (A 2 /8) exp (—0). 



There are field theoretical models that are able to describe statistical en- 
sembles of a discrete sets of elements and these models lead to the Liouville 
equation when a particular condition (called : self- duality) is fulfilled. We 
note that the Liouville equation can lead to the above current density equa- 
tion (Helmholtz-type) in an approximation where the function differs weakly 
from a background value. 

In the following we will discuss solutions of the soliton type of the Liouville 
equation A0 = (A 2 /8)exp (—(f)). It is commonly considered that the Liouville 
equation is solved by the formula given in terms of two arbitrary complex 
functions F and G 



(x,y) = - In 



F' (z) G (z*) 



[l + f 6 F(z)G(z*)\ 



where z = x + iy and * is the complex conjugate. However this form of the 
solution is too general to be useful. More popular is the particular form derived 
from this one, 

4> = In [cosh (kx) + e cos (ky)\ 

solution of A0 = (1 — e 2 ) k 2 exp (— 20). This describes the cat's eye vortex 
chains of the Kelvin-Stuart stream function, and the magnetic flux function in 
the chains of magnetic islands (Finn and Kaw). 

Our main objective is to identify a type of solution of the Liouville equation 
which is more localised on a magnetic surface compared with the profile of a 
magnetic island of the same helical symmetry. For this the factor k must be 
independent of the local q: the spatial extension of the density perturbation 
at a magnetic surface and the helical symmetry of that magnetic surface must 
be independent. This is clear in the case of the Liouville (and sm/i-Poisson) 
equation, since it is conformally invariant. We cannot be sure that the solution 
we are looking for can be of the form shown above. The correct approach is 
to find any solution of the Liouville equation. 

We now describe the systematic way of obtaining solutions of the Liouville 
equation on periodic domains, taking y as the poloidal and x as the radial 
coordinate. The method (Tracy et al.) consists of taking the Liouville equation 
as the limit of the sm/i-Poisson equation: 

A 2 = A 2 exp (-0) , V = <p-(3 

and let f3 — > oo. Then the sm/i-Poisson equation becomes the Liouville equa- 
tion after taking A — > A. The sm/i-Poisson equation is exactly integrable on 
periodic domains since it possesses a pair of Lax operators (Ting, Chen and 
Lee). The eigenvalue problem for the Lax operators identifies a spectrum of 
complex values where the two Bloch functions are not independent (for one 
eigenvalue we have only one eigenfunction). These nondegenerate eigenvalues 
are called main spectrum and represent branching points in the two-sheeted 
Riemann surface associated to the Wronskian of the two periodic solutions. It 
is shown that the problem of finding the unknown function <p is mapped on 
the problem of motion of a set of functions (auxiliary spectrum) on this Rie- 
mann surface. The equations of motion are nonlinear but they can be solved 



exactly, using the Abel transform. In this transformation the two-sheeted Rie- 
mann surface is mapped onto a compact hyperelliptic Riemann surface and the 
equations of motion are linearised as rotations along the cycles on this complex 
curve. The topology of the hyperelliptic Riemann surface is that of a sphere 
with a number of handles (the genus of the surface) which is given by the ge- 
ometry of the cuts in the complex plane needed to uniformize the first Riemann 
surface. Therefore the genus is actually determined by the main spectrum of 
the Lax linear operator, or, in concrete terms, by the boundary condition we 
require for the solution. The linear equations of motion are integrated leading 
to phases, which are linear combinations of the original variables, x and y. 
The number of phases is the number of type A cycles on the surface (like the 
short turn on a torus), or, the genus ; so, the number of phases is again de- 
termined by the main spectrum. We have to return to the original framework, 
and this represents the Jacobi inversion problem. It is solved exactly in terms 
of Riemann theta functions. The exact solution of the sin/i-Poisson equation 
in terms of Riemann theta function, 6, is 



<f)(x,y) = 2 In 



0(1+|1)' 



0(1) 

where 1 = k x x + k y y + 1 , 1 is a vector of constants, initial phases, and 



The physical content of the problem is in the square matrix C whose dimen- 
sion is half the number of eigenvalues in the main spectrum. The matrix C 
is obtained from integrals of a basis of differential one-forms defined on the 
hyperelliptic Riemann surface along the basis of closed paths (cycles). These 
integrals can be converted into integrals along closed paths on the plane of the 
spectral variable, around cuts or crossing these cuts. The geometrical aspect 
of this conversion is numerically complicated due to the jumps of the phases 
of the complex integrand at crossing the cuts. However, the symmetries of the 
main spectrum allows us to use general forms of the matrix 

C« = W N - 2j+1 C* N _ j+1 , j < N/2 

°jj — u iV-i+lj 

A particular choice of the entries of C (which obeys the symmetries) corre- 
sponds, physically, to a particular form of the boundary condition assumed for 
0, on a linear section of the periodic domain. 

The solution of the Liouville equation can be obtained from that for the 
sinh- Poisson equation in a certain limit. This limit has been translated into a 
particular distribution of the functions of the auxiliary spectrum (Tracy et al). 
For any (x, y) there are three classes according to the positions relative to the 
inversion circle, which is given by \E\ 2 = A 4 /256, E being the spectral variable. 



Solution of Liouville equation 




Figure 1: The streamfunction solution of the Liouville equation. 

First, one notes that the discrete points of the main spectrum are situated in 
certain positions around this circle: (1) there are N inversion pairs, (Ej, Epf+j), 
in the set E u E 2 , E 2N with: E N+j = A 4 / (25QE*) for j = 1, N. (2) there 
are M pairs E 2N +i, • -, E 2N+2M such that E 2N+k = X 2 a k /16, E 2N+M+k = 
(A 2 /16) /al, with a k independent of A. From each pair of the eigenvalues of 
the set (2) (situated near the inversion circle) there are M auxiliary functions, 
which scale as A 2 . The rest of the points of the auxiliary spectrum are devided 
into two classes. The first contains the auxiliary functions which are outside 
the inversion circle and are independent of A. The second are defined inside the 
inversion circle and are scaled as A 4 . In this way, A 2 exp (—(f)) is independent 
of A. This classifications cannot be directly employed, but they suggest a 
particular choice for the points of the main spectrum. These studies, mainly 
numerical, are still in progress. 

A conclusion can be drawn at this stage: the solution exhibits a localised 
perturbation on the poloidal direction, while the helical symmetry is still that 
given of the q of the surface. This solution is solitonic and therefore is sta- 
ble (we still have to clarify the effect of the limiting procedure on the set of 
invariants when passing from sm/i-Poisson to Liouville). The fact that this 
kind of solution is an attractor comes from the general property of the plasma 
state: the self-duality (leading to the Liouville equation) corresponds to the 
extremum of the action of the system of discrete elements. 
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